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Abstract 



We present a model improving the two-angle model for interphase chromatin 
(E2A model). This model takes into account the cylindrical shape of the 

histone octamers, the HI histones in front of the nucleosomes and the ver- 
tical distance d between the in and outgoing DNA strands. Factoring these 
chromatin features in, one gets essential changes in the chromatin phase di- 
agram: Not only the shape of the excluded- volume borderline changes but 
also the vertical distance d has a dramatic influence on the forbidden area. 
Furthermore, we examined the influence of HI defects on the properties of 
the chromatin fiber. Thus we present two possible strategies for chromatin 
compaction: The use of very dense states in the phase diagram in the gaps 
in the excluded volume borderline or missing HI histones which can lead to 
very compact fibers. The chromatin fiber might use both of these mecha- 
nisms to compact itself at least locally. Line densities computed within the 
model coincident with the experimental values. 
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Introduction 

The nucleosome is the basic repeat unit of chromatin m in all eucaryotic 
organisms. It consists of a cylindrical-shaped histone octamer and a stretch 
of DNA which is wrapped around the histone complex approximately 1.65 
times. The histone octamer consists of four pairs of core histones (H2A, 
H2B, H3 and H4) and is known up to atomistic resolution ([21 |3|. The 
nucleosomes are connected by blank DNA strands and together with these 
linkers they form the so-called 30nm fiber. The histone HI (and the variant 
histone H5 with similar structure and functions) is involved in the packing 
of the beads on a string structure into the 30nm chromatin structure. To 
do so it sits in front of the nucleosome keeping in place the DNA which is 
wrapped around the histone octamer and thus stabilizes the chromatin fiber. 
The nucleosome provides the lowest level of compaction and, furthermore, it 
is important in the regulation of transcription. Several enzymes can change 
the position of the nucleosome (4) to make the genetic information which is 
held within the nucleosome core particle accessible. 

The compaction of the DNA plays a very important role in modern bio- 
physics since it has a total length of some meters but has to fit into a 
nucleus of some microns. The degree of compaction depends on the salt 
concentration (|5]) and on the presence of linker histones (l6|). The presence 
of the linker histones leads to the formation of stem-like structures which 
are formed by the incoming and outcoming DNA string in front of the nu- 
cleosome. At low salt concentration a lOnm structure is formed which has 
the shape of beads on a string ([5$ whereas at high salt concentrations the 
chromatin fiber is much more compact and has a diameter of 30nm (JT)) . 

The chromatin structure is still not completely understood (jH [3 [9l [10]). 
There are competing models for its structure: zigzag ribbon models ([H 
[TT1 [T2l [T3l [HI, helical solenoid models ([5] [151 [T6|), superbeads or simply 
having no regular structure (fT7|l . A crystal structure of a tetranucleosome 
has been revealed (jlip and used to construct a model for the 30nm fiber 
which resembles a zigzag ribbon that twists or supercoils. The chromatin 
fiber has been investigated by electron cryo-microscopy ((6l [T8]) , atomic force 
microscopy (|19| I^U)) neutron scattering and scanning transmission electron 
microscopy (|2ip . Beyond the 30nm level chromatin is poorly understood. 

The two-angle model was introduced by Woodcock et al. (I12p to describe 
the geometry of the 30nm chromatin fiber. It has been shown that the 
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excluded volume of the histone complex plays a very important role for the 
stiffness of the chromatin fiber (22) and for the topological constraints during 
condensation/decondensation processes (|23p . In (|24p a rough approximation 
of the forbidden surface in the chromatin phase diagram was given. In 
a previous work of ours (f25]) we answered questions concerning the fine 
structure of the excluded volume borderline which separates the allowed 
and forbidden states in the phase diagram with the basic assumption of 
spherical nucleosomes and no vertical shift between in and outgoing strand. 
Here we present a Ramachandran-like diagram for chromatin fibers with 
cylindrical nucleosomes for a new extended model and furthermore discuss 
the influence of a vertical shift between the linkers due to HI histones and 
the volume exclusion of the DNA. 

We first present the basic notations for the formulation of the Extended 
Two- Angle model (E2A). Then we give an algorithm for the generation of 
chains within the model and present the resulting phase diagram and end- 
to-end distance as well as radius of gyration results. 

Theory 

Extended two-angle model 

We extend the two-angle model by introducing a parameter for the vertical 
distance between the DNA strands in front of the nucleosome. Furthermore, 
we take the cylindrical excluded volume of the nucleosomes into account as 
well as the HI histones which fix the DNA linkers in front of the nucleosome. 
The HI histones themselves are taken to be random variables to allow for 
possible missing HI histones. 

Basic notations 

We start out and fix some basic notations to use for the formulation of 
the model. The nucleosomes will be characterized by the centers Ni £ 
and the orientations pi G of the nucleosomes, with i = 0,...,N and 
IIPjII = 1. A'^ is the length of the fiber. The linkers between the centers of 
two nucleosomes will be denoted by 6i := Ni — Ni-i with i = 1, ...,N. The 
length ||6j|| of the linkers will be a further input parameter of the model 
(opposite of the direction hi G of the linkers). Furthermore, the entry- 
exit-angle Ui G [0,-7r] between two consecutive linkers is defined by Oj := 
<(— 6j, 6j_|_i) with i = 1, A^ — 1 and the rotational angle f3i G [0, vr] between 
two consecutive orientations is given by j3i := <{pi-i,pi) with i = 1, N . 



chromatin: phase diagram and iuEuence of HI defects 



4 



Moreover, hi represents the distance along the orientational axis Pi-i 
from Ni-i to Ni due to the spatial discrepancy between in and outgoing 
DNA strand, hi can be expressed by the vertical distances di which the DNA 
covers by wrapping up itself around the histone complexes: hi = ^{di-i + di) 
with i = 1, ...,N. 

Thus a single chromatin strand within the two-angle-model is character- 
ized by the following set of variables 

i{'^i}ie{l,..,N-l},{f^i}ie{l,..,N}Ahi}ie{l,..,N},{\\f^i\\}i€{l,..,N})- 

The general rotational matrix TZ around an axis d = (01,02,03)* (with 
||a|| = 1) by an angle 7 with respect to the right-hand rule will be used in 
the following. It is given by: 



a 



C0S7 ~{- a^{l — C0S7) aia2{l — C0S7) — a3sin7 aia^(l — C0S7) + a2sin7 

"JZ^ ^ I a2ai{l — C0S7) + a3sin7 C0S7 + 1^2(1 ~ C0S7) a20,^(l — C0S7) — aisin7 

a^ai(l — C0S7) — a2sin7 a^a2{l — C0S7) + aisin7 C0S7 + a|(l — C0S7) 



Definition of the two-angle model 

A chromatin fiber within the framework of the extended two-angle model 
has to fulfil the following conditions for all i G {1, A^}: 

■\ ^ f k u \ ^ c ^ {-bi,bi+i) 
i) <{-hi,hi+i) = ai<^ cos(aj) = „- 

\\0i\\ \\Oi+l\\ 

a) <{pi-i,pi) = I3i<^ cos(/3i) = {pi-i,pi) 
in) \\Ni - Aj_i|| = \\bi\\ 

iv) {pi-i,bi)pi-i = hi = ^{di-i + di). 

These are illustrated in Fig [T} 

The first condition adjusts the entry-exit angle of nucleosome i to the 
given parameter a,. The second condition does the same for the rotational 
angle due to the DNA twist from nucleosome i — 1 to nucleosome i. The 
third condition fixes the distance of the two consecutive nucleosomes i — 1 
and i and the last condition adjusts the vertical distance along the local 
chromatin axis between the nucleosomes i — 1 and i. 



Construction of the fiber 

The construction of the fiber can be done using an iterative process. A 
further part of the model is the presence of a HI histone which is assumed 
to be present with probability p. 
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The first nucleosome center and its orientation are arbitrary. We chose: 









: 


|, Po=| 











The following vectors fulfil the conditions of the two angle model for the 
second nucleosome location and its orientation: 

; / -VWblP-h'i 

N, = No + ^\\bi\\^-hl{ \+hipo=i 
and 

Pl=7^flK_l with a = (1,0,0)*. 

Now we can calculate Ni^i and Pi+i in dependence of Ni, A^j-i, Pi and 
Pi-i. With 

Vi := -bi + {pi, -bi)pi 

and 

< ■■= '^;:^Jbli-dli (^) + h^+m (1) 
one gets the location of nucleosome i + 1 by 

Ni+i = Ni + v'i. 

oq is the angle between the projections of and — onto an arbi- 
trary plane orthogonal to pi. We need to calculate the dependence of this 
projected entry-exit-angle ao on the actual entry-exit-angle a. 

Note that ao was used as entry-exit-angle in some other publications ()13| 
[M]) but in this work it denotes only the projection of the real entry-exit-angle 
a. 

Using the law of cosine one gets 

f = b^ + 62^1 - 2bA+icos{a) . (2) 

Now we will use an affine transformation T to a new coordinate system 
{x,y,z) — > {x',y',z') in order to get a second relation for I. We shift the 
origin to Ni and rotate our old coordinate system so that pi corresponds 
to the new z-axis. Furthermore, the new x-axis has to coincide with the 
projection of —bi onto any plane orthogonal to pi. 
Obviously, 

i' = \\h+vr = m+v'n\' 



chromatin: phase diagram and iuEuence of HI defects 



6 



with 



( ^bj-{p„-b.y 



{Pi 



and 



v' ^ v'/ 



cos^aojA/fe'+i 



"-i+l 



b? 



h^+i - cos(ao)AA?+i -/iLi 



H+l 



This leads to 



i+l 



+ h 



-2/ii+i(i)i,-6,)-2cos(ao)y6MS^^6^y6^^^^. (3) 



By comparing Eqj2] and Eqjs] one gets eventuahy 

6i6j+icos(a) - hi+i{pi, -bi) 



cos(ao) 



^bi,-hi,^b]-{p,,-b.y^ 



with the boundary condition 



ao > Oimin = acos 



-2hh 



(4) 



due to non- vanishing di and dj+i. The calculation of Ni^i is complete, since 
we now know the dependence of ao on a and therefore one can use EqJT] 
to determine iVj+i. But one still has to calculate the orientation Pi+i of 
nucleosome A^i+i- Due to the fixation of the in and outgoing DNA strand 
by the HI histones this orientation can be calculated by a rotation around 
the following normalized axis a: 



a :- 



-bi+i - {Pi, -bi+i)pi 



Pi-\-i then follows by a rotation of pi around this axis: 

Pi+i = nf+'p^. 
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Methods 

Chromatin phase diagram 

First we determine the influence of the cylindrical excluded volume of the 
nucleosomes and a non-vanishing vertical distance between in and outgoing 
DNA strand on the phase diagram of chromatin. Both of these parameters 
have been neglected so far (|24| |25|) . To do so we made simulations (I26p of 
regular chromatin fibers and checked whether they fulfil the excluded volume 
conditions or not. We were able to plot our results in a Ramachandran- 
like diagram (cf. Fig. [s]) and thus find out which states of the whole phase 
diagram are forbidden by excluded volume interactions and which are not. 
The fibers we simulated for this part were regular, i.e. all linker lengths, 
entry-exit-angles, rotational angles and hi were fixed for a certain strand. 

The cylinders had a height of 16.8(bp), and a diameter of 33.0(bp) ac- 
cording to ()27p . They were orientated by using the vectors pi above. More- 
over, we assumed a DNA diameter of 6.6(bp), a twist length of 10.2(bp), 
a mean linker length of 63 (bp) and that 1.65 turns of DNA are wrapped 
around the histone octameres. 

Fibers with HI defects 

Furthermore, we made Monte Carlo simulations of chromatin fibers with HI 
defects, i.e. some of he HI histones were missing. We used the two-angle- 
model with some fixed parameters (see above), which reflect the probable 
mean values within the cell. The only interaction potential is the hard core 
excluded volume. In this mean field-like approach we neglect (j28l) thermal 
fluctuations and thus assume to be above the 0-point to concentrate on the 
interaction between HI defects and volume exclusion. 

For a certain nucleosome iVj the defect probability p gives the chance of 
a missing HI histone. If the histone is missing the in and outgoing DNA 
strand are no longer fixed in front of the nucleosome but instead are arbitrary 
with respect to the excluded volume interactions of the chromatin strand 
(cf. Fig. [2]). Thus we get results for the mean end-to-end distance and the 
mean radius of gyration of fibers with various defect probabilities: p = 0.00, 
p = 0.01, p = 0.05, p = 0.10 and p = 0.30. For these simulations we fixed 
the entry-exit-angles Oi to 40 degrees, the rotational angles Pi to 36 degrees 
and hi = 7 (bp). 
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Results 

Phase diagram 

The colored lines in Fig. [3] represent the phase transition between allowed 
and forbidden states. All states below the corresponding line are forbidden, 
those above it are allowed. The states near the excluded volume borderline 
are the most interesting of the phase diagram since they are the most com- 
pact ones (cf. Fig. |4]). Therefore, the gaps in the borderline might be used 
by the fiber to become (at least locally) very dense. 

There is another borderline at the left side of the diagram which prevents 
a from getting smaller than some minimal value amin{h) which depends on 
h. This arrow-like structure can be seen best in Fig. |6] It shifts towards 
larger amin{h) with increasing d. The gap in this line is a further conse- 
quence of the cylindrical excluded volume and can not be seen in the phase 
diagram for spherical nucleosomes (|25p . 

As a consequence of the cylindrical instead of spherical excluded volume 
of the nucleosomes the shape of the peaks in the phase diagram is changed: 
Their top shows a wedge- like shape due to the edges of the cylinders. With 
increasing entry-exit-angle a there is more space between the nucleosomes 
which leads to a larger variety of allowed rotational angles /3 and thus to the 
missing tip at the top of the peaks. This effect gets weaker with increasing 
d: The edges which cut the peaks become more parallel to the a- axis. 

With increasing vertical distance between in and outgoing DNA strand 
there is also more space between consecutive nucleosomes which leads to 
a decrease of the borderline. More and more states become accessible and 
with h = 5. 5 (bp) the borderline almost vanishes. The natural mean value 
of h is approximately h = 2.8nm = 8.4(bp) due to 1.65 turns of DNA with 
a diameter of 2.2nm. Lower h values might occur where the DNA has less 
turns around the histone complex. 

We furthermore examined the line density (cf. Fig. [5]) and the radius of 
gyration (cf. Fig. |6]) of regular chromatin fibers (length 500 nucleosomes) 
along the phase diagram. The most compact states can be found near the 
excluded volume borderline. The line densities we found in our simulations 
coincide with experimental values (j29D. Increasing d decreases the line den- 
sity and increases the radius of gyration. 

HI defects 

We also investigated the influence of missing HI histones on the mean 
squared end-to-end distance (Fig. [8| and the mean squared radius of gy- 
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ration (Fig. |9]) of chromatin fibers. The parameter p gives the defect prob- 
abihty in this section. One can clearly see that even very small defect rates 
of some percent have a huge effect on the compaction of chromatin: Both 
mean squared radius of gyration and the mean end-to-end distance decrease 
rapidly if one allows only a few HI defects. Without HI defects {p = 0) we 
get an ideal chromatin fiber within the restriction of the extended two-angle 
model. This ideal fiber reflects the properties of the 30nm strand only on 
small length scales. Therefore, the increase of the compaction due to defects 
will probably be not as strong as implicated by our results. Nevertheless 
missing HI histones might contribute to chromatin compaction and DNA 
accessibility for transcription purposes at the same time since one can see 
from Fig. [7] that although the fiber gets compacter some very open parts 
appear. One can also see here the increasing disorder with increasing p. 

Discussion 

The compaction of chromatin is still an open question. A polymer of a total 
length of two meters has to fit into a tiny cell nucleus of some microns. We 
showed two possible strategies for the fiber to deal with this task. The fiber 
might use gaps in the phase diagram, i.e. very dense states to compact parts 
of itself. Furthermore, we showed that missing HI histones might supply a 
further contribution to the compaction of the fiber. These HI defects might 
play a very crucial role in the task of chromatin compaction and at the same 
time serve the transcription of the DNA by opening locally the chromatin 
fiber. 

Moreover, we developed the ordinary two-angle model further to our E2A 
model, which is much more detailed and thus appropriate to model chro- 
matin at the 30nm level. 
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Figure 1: The figure shows the basic parameters of the E2A model: The 
entry-exit-angle Oj, the rotational angle Pi, the linker length bi and the 
vertical distance di between in and outgoing linker. We chose a large entry- 
exit-angle here to make the visualization clear. 




Figure 2: Shown is an example of a HI defect within the chromatin fiber. 
The upper strand and the strand below the defect are regular. 
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Phase diagram of the chromatin fiber 




a [deg] 



Figure 3: Cut-out of the chromatin phase diagram (for different d). The 
states below the corresponding lines are forbidden due to excluded volume 
interactions. With increasing d more and more states become accessible to 
the fiber. 
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End-to-end-distance [nm], d=Onm 




Figure 4: The end-to-end distance of regular chromatin fibers along a cut- 
out of the phase diagram with fixed fiber length (A^=500 segments) and 
fixed (i=0.0(bp). The solid black line represents the corresponding phase 
transition. 
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Line density [#nucleosomes/1 1nm], d=Onm 




a [deg] 

Line density [#nucleosomes/1 1 nm], d=2nm 
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Line density [#nucleosomes/1 1 nm], d=1 nm 




a [deg] 

Line density [#nucleosomes/11nm], d=3nm 




a [deg] 



Figure 5: The line density of the nucleosomes within a regular chromatin 
fiber is large for small However the comparison with the phase dia- 
gram shows, that some of these states are forbidden due to excluded volume 
interactions. 
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Radius of gyration [nm], d=O.Onm Radius of gyration [nm], d=1 .Onm 




a [deg] a [deg] 



Radius of gyration [nm], d=2.0nm Radius of gyration [nm], d=3.0nm 




a [deg] a [deg] 



Figure 6: The radius of gyration of regular chromatin fibers along a cut- 
out of the phase diagram with fixed fiber length of A^=500 segments. The 
compaction of the fibers decreases strongly with increasing d. 
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Figure 7: Chromatin fibers with different defect probabihties p. At p = 0.30 
the regular structure of the 30nm strand is almost completely collapsed. 
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End-to-end distance for chromatin fibers with H1 defects 




genomic distance [kbp] 



Figure 8: The mean end-to-end distance for chromatin fibers with HI de- 
fects. With increasing defect probabihty p the length of the fibers decreases 
rapidly. HI defects might play a crucial role for chromatin compaction. 
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Squared radius of gyration for chromatin fibers with H1 defects 




5 10 15 20 25 30 35 40 45 50 55 60 65 
genomic distance [kbp] 



Figure 9: The squared radius of gyration for chromatin fibers with differ- 
ent defect probabihties p. With increasing number of HI defects the fiber 
becomes much more compact which could be an important mechanism to 
compact the chromatin fiber. 



